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Abstract. We present 4 very low mass stars radii measured with the VLTI using the 2.2pm VINCI test instrument. The 
observations were carried out during the commissioning of the 104-meter-baseline with two 8-meter-telescopes. We measure 
angular diameters of 0.7-1.5 mas with accuracies of 0.04-0.11 mas, and for spectral type ranging from MOV to M5.5V. We 
determine an empirical mass-radius relation for M dwarfs based on all available radius measurements. The observed relation 
agrees well with theoretical models at the present accuracy level, with possible discrepancy around 0.5-0.8 Mq that needs to be 
confirmed. In the near future, dozens of M dwarfs radii will be measured with 0.1-1% accuracy, with the VLTI, thanks to the 
improvements expected from the near infrared instrument AMBER. This will bring strong observational constraints on both 
atmosphere and interior physics. 
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1. Introduction 


Mass, radius, luminosity and chemical composition are the ba¬ 
sic physical properties of a star. For a given mass and chem¬ 
ical composition, theory can predict most of the stellar phys¬ 
ical parameters at a given age. Accurate stellar mass, radius 
and luminosity measurements thus provide a crucial test of 
our understanding of stellar physics. To be relevant, these pa¬ 
rameters muMte determined with an accuracy of the order 
of 1% ( Andersen 1991 ).The previously rather noisy Mass- 
Luminosity relation for M dwarfs has recently been greatly im¬ 
proved (Delfosse et al. 200C; Segransan et al. 2003[ ). By con¬ 
trast, the empirical Mass-Radius relation remains poorly con¬ 
strained, since it is based on the observations of just two of the 
three known M-type eclipsing binaries, plus three M dwarfs 
radii measured at the Palomar Testbed Interferometer (PTI) 


with accuracies of 2-4%, (Lane et al. 2001). 

In this paper we present direct angular diameter mea¬ 
surements of four M dwarfs with spectral types rang¬ 
ing from MOV to M5.5V obtained with the Very Large 
Telescope Interferometer (VLTI) and two 8-meter-telescopes 
( [Glindemann et ah 2001).The first two sections describe the 
observations and the data reduction. The third section presents 
the visibility calibration and the angular diameter determina- 
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tion. In the last section we compare presently known M dwarfs 
radiis and masses to theoretical models. 


2. Observations 

G1887, G1205 and G1191 were observed during the first com¬ 
missioning run of the VLTI with 8-meter-telescopes in early 
November 2001, on the second and third nights after the “first 
fringes” with these telescopes. The observations of G1551 were 
obtained in February 2002. 

VINC I is th e commissioning instrument of the VLTI (Kervella 
et al. 2000). It uses single-mode optical fibers to recombine 
the light from two telescopes of the Paranal Observatory, and 
modulates the optical path difference around its zero position to 
produce interference fringes. This recombination scheme was 


first used in the FLUOR instrument ( Loude Du Foresto et al 


1998) and produces high precision visibility values, thanks to 


the efficient spatial filtering of the incoming beams and to the 
photometric monitoring of the filtered wavefronts. 

Transfer function calibrators (Table were selected based on 
their spectral type, distance to the target and known diameter 
(Cohen et al. 1999). Calibration tests show that a +1 cr shift 


in the calibrator angular diameter values modify the measured 
M dwarf radii by less than 0.8 cr. 

The M dwarfs science targets and their calibrators were 
observed in series of a few hundred interferograms, recorded 
at a frequency of 500-700 Hz. Due to many down times 
associated with the technical commissioning of the array, the 
observations of calibrator stars were not possible to the extent 
usually required to monitor extensively the transfer function. 
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Table 1. Calibrators and instrumental+atmospherical transfer 
function. Calibrator uniform disk angular diameters come from 
Cohen et al. (1999). Internal error and external error on the 
transfer function are given for each calibrator. Fomalhaut an¬ 
gular diameter was redetermined using the VLTI. (*) calibrator 
was not used due low visibilities or poor accuracy on their di¬ 
ameters 


Targets 

Calib. 

Spec. 

Type 

Diameter 

0(/o[mas] 

Transfer function 

l^a-Fil ^ tot iP " stat ! ^ cal ) 

GJ205 

39-Eri 

K3III 

1.81±0.02 

0.478 

0.005 

(.002/.005) 


HD36167* 

K5III 

3.55±0.06 

0.486 

0.037 

(.006/.036) 

GJ887 

HR8685 

MOIII 

2.01±0.02 

0.448 

0.005 

(.003/.004) 


Fomalhaut* 

A3V 

2.13±0.06 

0.433 

0.013 

(.004/.013) 

GJ191 

39-Eri 

K3III 

1.81±0.02 

0.420 

0.005 

(.004/.004) 


Gam02 Vol* 

KOIII 

2.45±0.06 

0.411 

0.014 

(.006/.013) 

GJ551 

HDl 10458 

KOIII 

1.66±0.02 

0.509 

0.006 

(.004/.004) 


Nevertheless, the calibration measurements proved sufficient 
to estimate the transfer function with a good accuracy. The 
global interferometric efficiency of the VLTI mostly depends 
on two parameters; the atmospheric coherence time, which 
varies over timescales sometimes as short as a few minutes, 
and the instrumental efficiency which changes over timescales 
of a week. Our effective calibration rate oversamples the 
instrumental efficiency variations by orders of magnitude, 
but risks undersampling the atmospheric coherence time. 
Fortunately, the Paranal observatory site is heavily monitored, 
and we could rely on cont inuous measurement o f the atmo¬ 
spherical coherence time ( Sandrock et al. 2000 ) to discard 
data intervals when the coherence time changed significantly 
between the observation of the calibrators and science targets. 


3. From fringes to visibilities 


We used a customised version of the VLTI/VINCI data reduc¬ 
tion pipeline, whose g eneral principle is based on the origi¬ 
nal FLUOR algorithm ( Coude Du Foresto et al. 1997 ). Despite 
the high modulation frequency of the fringes, many recorded 
interferograms present a differential piston signature, as well 
as strong photometric fluctuations. The latter are expected for 
large apertures without adaptive optics correction, as individual 
speckle come in and out of the fiber input aperture. To over¬ 
come the noise amplification caused by low photometry data 
points, which could strongly bias the resulting visibilities, in¬ 
terferograms were calibrated by the average value of the pho¬ 
tometry over the fringe packet. The two calibrated output inter¬ 
ferograms were then subtracted to remove residual photometric 
fluctuations (Kervella, in prep). 

In parallel to the FLUOR algorithm, based on Fourier anal¬ 
ysis, we implemented a time-frequency analysis (Segransan 
et al. 1999) based on the continuous wavelet transform (Farge 


1992). Instead of the projection onto a sine wave of the Fourier 


transform, the wavelet transform decomposes it onto a func¬ 
tion, ie. the Wavelet, that is localised in both time and fre¬ 
quency. We use as a basis the Morlet Wavelet, a gaussian en- 
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Fig. 1. Wavelet Spectral density of a pistonned interferogram. 
The effect of strong differential piston is to distort the interfer¬ 
ogram fringe peak in the time-frequency space. Part of the in- 
terferogram’s energy is spread out in both time and frequency, 
preventing a good measurement of the total energy. 


veloppe multiplied by sine wave. With the proper choice of the 
number of oscillations inside the gaussian envelope, the Morlet 
wavelet closely matches a VINCI interferogram. It is therefore 
extremely efficient at localizing the signal in both time and fre¬ 
quency. 

As illustrated on figure ^ differential piston strongly affects the 
amplitude and the shape of the fringe peak in the wavelet power 
spectrum. We therefore select on the shape of fringe peak in the 
time-frequency domain to efficiently reject “pistonned” inter¬ 
ferograms. We then derive visibilities from the wavelet spectral 
density, after removing the residual photon and detector noise. 


4. Angular Diameter & Limb Darkening 

Stellar angular diameters are usually quoted for both a uniform 
disk and a limb-darkened disk, though only the latter is directly 
comparable to actual stellar mod els. The monochr omatic visi¬ 
bility of a limb darkened disk is (|Davis et'ffi^ 2000): 


Va = 


dfilin)nJoi^nBOLDl ^ j 

fo 


( 1 ) 


where / (/j) is the surface brightness, Jo is the zero'^' order 
Bessel function, B the projected baseline, A the wavelength, 
and Gld the limb darkened diameter (in radian). 

The VINCI instrument has no spectral resolution, and a 
bandpass defined hy a K filter (2-2.3/im). For marginally re¬ 
solved disks {BOjA « 1), Eq. |] remains approximately valid 
for polychromatic measurements, with A then being the ef¬ 
fective wavelength of the system. It differs from the mean 
wavelength of the filter through the wavelength-dependency of 
fiber coupling, combiner transmission, and the stellar spectrum. 
Since interferograms are Fourier transforms of the instrumental 
bandpass, they provide, in principle, an accurate measurement 
of the effective wavelength, obtained as the average barycen- 
ter of the fringe peak in the spectral density space. We com¬ 
puted such effective wavelengths for each star and compared 
them to the one given by the VINCI instrument modelisation 
{Aeff - 2.195 + 0.010/im). Our estimations of Aeff present a 
large scatter which is due to atmosperical differential piston and 
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Fig. 2. Measured visibilities (data points) and best-fit model 
(curves, equation ^ for GJ 205 (plain), GJ887 (dots), G1191 
(short dash), and G1191 (long dash) 


photometric fluctuations. We thus used the VLTl/VINCI effec¬ 
tive wavelength provided by ESO and the effective temperature 
of each star to compute an accurate Agff for each target and cal¬ 
ibrator. 

Due to their relatively short baselines, and limited uv- 
coverage, these observations cannot discriminate between an 
unphysical uniform star and a more realistic slightly larger 
limb-darkened star. To correct for limb-darkening, we there¬ 
fore have to use a theoretical limb-darkening law. We adopt 


the non-linear limb-darkening coefficients from Claret (2000), 
Eq. based on atmospheric models from Allard et al. (1997). 


We us e gravities and effective temperatures from Baraffe et al. 
( 1998 ) isoc hrones at 5 G yr (Table ||) to select the apropriate 
entry in the Claret (200C) tables. 



( 2 ) 


In Equation the at are the limb-darkening coefficients and 
jj. - cos{y), with y the angle between the line of sight and the 
emergent ray. The limb-darkening correction that we apply 
is in principle model-dependent. Eortunately, it is also quite 
small for M dwarfs, only 1-2% (Table The correction itself 
is at best marginally relevant at our current accuracy level 
for individual diameters , and its theoretical uncertainties can 
for now be safely neglected. Euture observations with longer 
baselines, at shorter wavelengths, and with a more accurate 
calibration, will need to consider the issue more carefully, but 
they will also directly measure the limb-darkening law for the 
closest M dwarfs. 


5. Mass-Radius relation 

An accurate empirical mass-radius relation is an essential con¬ 
straint on stellar interior structure, evolutionnary models and 
atmospheric physics. The interior structure is largely deter¬ 
mined by the equation of state, whose derivation for very low 


mass stars, brown dwarfs, and planets involves the complex 
physics of strongly correlated and partially degenerated quan¬ 
tum plasma (Chabrier & Baraffe 200C). 

To obtain the empirical mass-radius relation we need to con¬ 
vert angular diameter into linear radius, which is easily done 
thanks to accurate HIPPARCOS parallaxes for all sources, and 
we need accurate masses. Eor these single stars masses can only 
be estimated, from IR photometric measurements and accu¬ 
rate mass-luminosity relations. Eortunately, the K-band mass- 
luminosity relation has very little intrinsic dispersion for M 
dwarfs . We estimated the masses of the stars listed in Table 
fusing an update (Segransan et al. 2003) of the (Delfosse et al 


2000) empirical K band mass-luminosity relation, now based 
on 27 accurate masses and luminosities. The data dispersion 
around the average empirical relation corresponds to a mass 
error of ~ 5%. Masses and radii are summarized in Table 

Eigure ^ compares the empirical r adii & masses with 5 Gyr 
and 1 Gyr theoretical isochrones from Baraffe et al. ( 1998| ). The 
Baraffe et al.| ( |1998 ) models reproduce the observation fairly 
well in the sampled range, between 0.65 and 0.12 Mq. At a 
more detailed level though, one notices that the models under¬ 
estimate the radii for YY Gem, V818 Tau and GJ 205, with 
masses in the 0.5-0.65Mo range. The discrepancy is highly sig¬ 
nificant for the eclipsing binaries, as extensively discussed by 
Torres & Ribas (2002) and is more marginal for long base¬ 
line interferometry data. There is an indication in Eig. H that 
the model reproduces the observations well below 0.5 Mq, and 
only become discrepant above that value. If real, this suggests 
that the shortcomings of current models have to be searched in 
the energy transport (convection description, opacities), rather 
than in the equation of state (EOS). As their mass decreases, 
stars at the bottom of the main sequence have increasingly 
simple transport properties (they are fully convective below 
~0.3 Mq), and an increasingly correlated and degenerate EOS. 
EOS shortcomings are therefore expected to show up most 
prominently at the lowest masses, and transport problems at 
higher masses. That result obviously needs confirmation from 
additional data points and from more accurate measurements, 
but the consistency between the eclipsing binary and the re¬ 
solved single stars is comforting, as these two datasets were 
obtained through completely independent methods. 

In addition to providing physical radii, the angular radii can be 
combined with integrated flux measurement to obtain the ef¬ 
fective temperature from the Stefan-Boltzmann law. Up to now 
that process has typically been inversed, to derive radii for stars 
whose effective temperatures were obtained from comparison 
with model spectra (e.g. Leggett et al. (2000)). Amongst the 
sources in Table ^ only GJ 699 has an integrated flux mea¬ 
surement (Leggett et al. 2000), and we thus have to rely on the 
bolometric correction polynomial fit of (Leggett et al, 2000) to 
estimate it for the other targets. 

In spite of these increased uncertainties on the bolometric 
flux, the resulting effective temperatures have fairly small error 
bars (Table With increasingly precise direct angular radii, 
and with integrated flux measurements for all targets, this test 
can become a very discriminating validation of model atmo¬ 
spheres, with potential accuracies at the ~ 1% level on both 
Tg/y and log(g). The present diameter measurements, which 
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Table 2. Measured angular diameter for the target stars. The uniform disk diameters for GJ 699, GJ 15A, GJ 411, GJ 380, 
GJ 105A are from Lane et al. (2001). For consistency with our measurements, we recomputed their limb darkened diameters, 
masses and radii. The masses for G1 380 and G1105A are based on the Baraffe et al. (1998) M-L relation, with an arbitrary 5% 
error bar. 



Spectral 

Photometry & Masses 

Limb Darkening, K band (models) 

Diameter [mas] 

Radius [Rq] 

Atm. prop. 

Object 

Type 

Mk 

M/Mo 

g 

Teff 


^3 d,^ 

SuD 

Sld 

o-g 

R 

o-R 

Teff 


g 

0's 

GJ205 

M1.5V 

5.09 

0.631±0.031 

4.70 

3894 

1.11 -1.11 

0.92 -0.31 

1.124 

1.149 

0.11 

0.702 

0.063 

3520 

170 

4.54 

0.06 

GJ887 

M0.5V 

5.79 

0.503±0.025 

4.80 

3645 

1.61 -2.35 

2.00 -0.68 

1.366 

1.388 

0.04 

0.491 

0.014 

3626 

56 

4.76 

0.03 

GJ191 

MIV 

7.08 

0.281±0.014 

4.98 

3419 

1.76 -2.72 

2.39 -0.82 

0.681 

0.692 

0.06 

0.291 

0.025 

3570 

156 

4.96 

0.13 

GJ551 

M5.5V 

8.80 

0.123±0.006 

5.19 

3006 

1.94 -2.80 

2.39 -0.81 

1.023 

1.044 

0.08 

0.145 

0.011 

3042 

117 

5.20 

0.23 

GJ699 

M4V 

8.21 

0.158±.008 

5.11 

3193 

1.87 -2.88 

2.54 -0.88 

0.987 

1.004 

0.04 

0.196 

0.008 

3163 

65 

5.05 

0.09 

GJ15A 

M2V 

6.27 

0.414±.021 

4.87 

3541 

1.66 -2.48 

2.14 -0.73 

0.984 

1.000 

0.05 

0.383 

0.02 

3698 

95 

4.89 

0.06 

GJ411 

M1.5V 

6.33 

0.403±.020 

4.88 

3533 

1.67 -2.50 

2.16 -0.73 

1.413 

1.436 

0.03 

0.393 

0.008 

3570 

42 

4.85 

0.03 

GJ380 

K7V 

4.77 

0.670±.033* 

4.65 

4106 

1.09 -1.01 

0.83 -0.28 

1.268 

1.155 

0.04 

0.605 

0.02 

- 

- 

4.70 

0.03 

GJ105A 

K3V 

4.17 

0.790±.039* 

4.56 

4603 

0.86 -0.53 

0.38 -0.13 

0.914 

0.936 

0.07 

0.708 

0.05 

- 

- 

4.63 

0.05 


have similar precisions to those obtained at PTI (Lane et al. 


2001 ), have been obtained during the commissioning phase of 
the VLTI interferometer, with a recombiner that was specified 
to validate and debug the interferometer, and not for an opti¬ 
mum scientific output. With the improvements e xpected from 
the near infrared scientific recombiner, AMBER (|Petrov et al. 


2001 ), from adaptive optics on the unit telescopes, the VLTI 
will measure dozens of M dwarfs with 0.1-1 % accuracy. This 
will bring very strong observational constraints on models of 
both stellar interiors and atmospheres. 
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